Nuclear response for the Skyrme effective interaction with zero-range tensor terms. 

II. Sum rules and instabilities. 
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The formahsm of hnear response theory for Skyrme forces including tensor terms presented in 
article [J] is generalized for the case of a Skyrme energy density functional in infinite matter. We 
also present analytical results for the odd-power sum rules, with particular attention to the inverse 
energy weighted sum rule, M_i, as a tool to detect instabilities in Skyrme functionals. 
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I. INTRODUCTION 

The Energy Density Functional (EDF) method is a tool 
of choice to perform systematic calculations of binding 
energies and one-body observables in the region of the 
nuclear chart that ranges from medium to heavy mass 
atomic nuclei from drip line to drip line 0. This ef- 
fective approach relies on a limited number of universal 
parameters, usually fitted on experimental data (observ- 
ables) Q along with properties of infinite nuclear mat- 
ter (pseudo-observables) extracted from experimental re- 
sults or derived from realistic models @. 

In its general formulation, the EDF is the sum of dif- 
ferent terms that depend only on products on one-body 
densities weighted by coupling constants, which in gen- 
eral can also depend on the local densities themselves. 
Although several functionals exist on the market 
the most often used is the one derived from the effective 
Skyrme interaction Building all possible combina- 
tions up to quadratic terms in densities together with 
the conservation of some general symmetries (see ref. @ 
for a detailed discussion), one obtains 28 free coupling 
constants that can be reduced to 14 imposing that 
the functional is derived from an effective force. This 
requirement is not only adopted to simplify the opti- 
mization procedure used to determine the values of the 
coupling constants, but it is also mandatory in order to 
use methods which go beyond the mean-field to avoid 
self- interactions and self-pairing [TT| - [l3j . Due to this 
additional requirement, it turns out that the standard 
Skyrme EDF usually adopted in mean-field calculation 
is not flexible enough to improve the level of accuracy 
in describing set of available experimental data. For this 
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reason, other terms are now considered as 3-body |14l | 
and tensor forces, for example. 

The determination of accurate values for the coupling 
constants of the Skyrme functional, even in its simplest 
form, it is quite an elaborated process since a good 
predictive power is only possible from a pertinent and 
well-chosen set of observables or pseudo-observables. Al- 
though this is still possible for the time-even part of the 
functional Q , there is not yet a consensus on how to pro- 
ceed in order to constrain the time-odd terms. Actually, 
it is not clear at all which observables could be used for 
this purpose, and the corresponding parameters are not 
explicitly constrained, but merely indirectly determined 
by the time-even part through simple mathematical rela- 
tions. For this reason, a vast area can be explored in the 
parameter space and possibly ends up in some region of 
instabilities as discussed further in this article. 

One of the first method used to fix some of these 
terms can be found in the work of Van Giai and Sagawa 
p^ . where they adjust the Landau parameters on val- 
ues obtained from realistic forces. Furthermore, from 
the theory of Landau-Migdal for quasi-particles, one 
can derive set of sum-rules for Landau parameters [l6j 
that should be fulfilled, otherwise the system could pass 
through different phase transitions according to the dif- 
ferent spin/isospin channels (for instance ferromagnetic 
instabilities in spin channels). The Landau-Migdal ap- 
proach is valid for quasi-particles interacting near the 
Fermi surface with transfer momentum that goes to zero 
- a situation that corresponds to the so-called long wave- 
length limit. Thus it is not able to predict possible 
instabilities that occur at non-zero transferred momen- 
tum q, with the appearance of domains with typical size 
A ~ ^Tr/q. The first example of such kind of instabil- 
ity was encountered and examined in details in an ar- 
ticle devoted to the study of effective mass splitting by 
Lesinski et al. [l3l , in the scalar-isoscalar channel of the 
SkP functional. It has been shown that when performing 
high-accuracy Hartree-Fock calculations (HF) of doubly- 
magic nuclei, the system converges towards an unphysi- 
cal configuration where protons are separated from neu- 
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trons. This observation has also been confirmed by RPA 
calculations in finite nuclei [isj . Another recent exam- 
ple of instability was found by Hellemans et al. [l^ in 
the vector channel of several Skyrme functionals. They 
have performed cranked-HFB calculations in ^^"'Hg and 
shown that for particular values of the time odd coupling 
constants the system can spontaneously polarize. 

To improve the existing functionals, it is therefore 
mandatory to find a tool which is able to detect these in- 
stabilities in all scalar (vector) /isoscalar (isovector) chan- 
nels. It has already been demonstrated by Lesinski et 
al. that the linear response formalism (LR) applied 
to the Skyrme energy functional could be used to predict 
the appearance of some finite-size instabilities in nuclei. 
However, only the central part of the Skyrme interac- 
tion was taken into account for the building of the lin- 
ear response. The same LR formalism for a case of a 
Skyrme interaction including tensor and spin-orbit term 
was studied by Davesne et al. [H, hereafter denoted as 
article I. In the present article, we extend the formalism 
of article I by expressing our results in terms of cou- 
pling constant of a general Skyrme functional. The main 
goal is actually to investigate the role of odd-power sum 
rules [2^ and show that they can be used as a simple 
and very powerful tool to detect instabilities in infinite 
symmetric matter. In particular we give the explicit ex- 
pression of the inverse sum rule M_i and demonstrate 
that a pole in the response function corresponds to a 
zero in the denominator of this inverse sum rules. This 
greatly simplifies the process of poles detection since we 
just have to find the roots of a real function. A more 
detailed analysis concerning the correspondence between 
finite size instabilities in finite nuclei and infinite matter 
will be the subject of a forthcoming article [2ll |. 

This work is organized as follows: in section|II]we sum- 
marize the different components of the EDF and recall 
the main steps of the LR formalism in nuclear matter 
presented in article I [H. In particular we give explicit 
expressions of the RPA responses in terms of the coupling 
constants of the Skyrme functional. We also establish 
the expressions of the first odd moments of the strength 
function in each channel. In section III, we present the 
results concerning the detection of instabilities and for 
completeness we also show the resulting Landau param- 
eters. Further possible developments are discussed in the 
conclusion. 



II. LINEAR RESPONSE 

A. Response functions and energy density 
functional 

The response functions, x'^'*^'^('^j Ct)i we are inter- 
ested in are formally defined as the response of the in- 



finite medium to external probes of the type Q^^'^'^-* = 
e''i '"j 0^^'^^'^^ where S (M) is the spin (its projec- 
tion along the z-axis), / the isospin and the operators 
0js,Mj) given in Table I. Following the notation of 
Garcia- Recio et al. [2^ . we have 
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where w, q are respectively the transferred energy and 
momentum, the sum is on all excited states \n) with en- 
ergy En and y is a quantification volume (see ref. [l^l 
for a detailed discussion) . Without any residual particle- 
hole, ph, interaction, the above expression reduces to the 
usual Lindhard function. Switching the ph interaction 
on, the response functions can be determined with the 
use of the RPA formalism. Such calculations have al- 
ready been published in the literature with Skyrme point- 
like interactions which incorporate only the central part 
[TtI . [2^ or the spin-orbit contribution [23] as well. More 
recently, in article I, we generalized the previous calcu- 
lations by taking into account the tensor part, which re- 
veals to be very important quantitatively. However, in 
article I we expressed the response functions with respect 
to the usual coupling constants of the Skyrme effective in- 
teraction : {ti,Xi},i = 0, ..,3, Wq and te,^o respectively 
for the central, the spin-orbit and the tensor part. In 
the present article we write them with an energy density 
functional (EDF) as a starting point. This has the great 
advantage to be more general in the sense that all the 
coefficients can be now considered as independent one 
from each other. The parameter space is thus enlarged, 
allowing by instance more fiexibility for the description 
of nuclei. In the context of forthcoming articles on insta- 
bilities, it will allow us to study precisely the role of each 
of these coefficients. To be specific, we consider hereafter 
the following EDF (see article I for notations) 



3 



-'Skyrme 



+C^-'{ptV ■Jt + Sf'^ xit) 



(2) 
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When the EDF is derived from a Skyrme interaction, 
the coupUng constants can be re-expressed, following the 
notation of article I, in terms of At and Bt coefficients. 
The coupling constant written At depend on the cen- 
tral and spin-orbit part of the interaction (i.e. Cf — A^ , 
C^" ^ Af", ^ Ai and Ct = Af, C^' = Aj-') 
and the ones written Bt depend on the tensor part (i.e 
C'^'^ = BY'^ and Ct = Bf), but can also contribute to 
the central part of the interaction {i.e. Cj = Aj + Bj 
and C(^* = A^^ +B^^). The expressions of the coupling 
constants as functions of the parameters of the interac- 
tion can be found in article I. The procedure used to 
obtain the residual interaction is then no longer based, 
as in article I, on the determination of the matrix ele- 
ments of the particle-hole interaction from the Skyrme 
one, but on the direct double derivation with respect to 



the one-body density of the EDF. The results concern- 
ing the residual interaction coming from the tensor part 
arc summarized in Appendix \^ while the response func- 
tions for infinite nuclear matter are explicitly written in 
Appendix |B] 

For completeness we also give the Landau parameters. 
To obtain their expression, we have to take the limit q — ^ 
and qi,2 -> 



Landau , 
ph ' 



kF,kF)= lim Vph{(l,(ii,c\2). (3) 



q— >0,qi,2— >kF 



According to ref . [ly, |25| , the most general form of the 
residual interaction in the Landau limit is 



V ^Landau 
ph 



5{vi - V2 



Hi Sab +%-H[ SabTa O f | (cOS 0)(4) 



— 1 h 

where Nr. = „ J, is the usual normalization fac- 
tor given here for the symmetric infinite nuclear matter, 
ki2 = (ki -k2) and Sab = 3(qi2 ■ ffa) (qi2 • cTb) - ctq • erf, , 
where the symbol qi2 indicates a vector of unitary length. 
One can express the product of momentum and Pauli 



matrices as (ki2 • cra)(ki2 • crt) = ^Sab + \o-a ■ ""b- It is 
important to notice that the H coefficients are functions 
of cos 9 and are only related to the tensor part of the in- 
teraction [l^, [26l - [29| . In our case the residual interaction 
in the Landau limit reads 



T/Landau/i^ \ 

Vph (k_F,kFj 



(0,0) 
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(0,0) 
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(0,1 



2,L TaOTb 



a ■ CTb 
(Tb + 
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(1,1), 



CTbTaO n 



(5) 



(Tb 



[2kj, - 2kj. cos e] 



[/Cj^Cq + kpCi f a O f b ] T2~ "^ai 



where the wj'^'^^ with j = 1,2 coefficients are given force we get the same values as in ref. [l6j . 
in Appendix [C] We checked that in the case of a Skyrme 
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B. Sum rules and moments of the strength function 

It should be noted that the quantity of interest is not 
directly the response function itself discussed in the pre- 
vious paragraph, but merely S'(")(q, w), usually called 
dynamical structure function, which is, at zero tempera- 
ture, proportional to the imaginary part of the response 
function at positive energies 



(6) 



On the other side, the k-moments, which are defined as 
moments per particle in infinite matter, read 



(7) 



is an analytic expansion (sec the next paragraph for ex- 
plicit expression). Both should coincide with very high 
accuracy. When it is not the case, this means that the 
dispersion relation is no longer valid or, in other words, 
that a pole occurs. Thus, a discrepancy between the dif- 
ferent expressions for the sum rules will be interpreted as 
the presence of a pole. 

We shall now enter the details and discuss separately 
the three important sum rules Mi,M3,M_i. Since the 
other cases can be obtained by switching off the appropri- 
ate coupling constants, the general case with the tensor 
will be considered only. As stated previously, all the ex- 
pressions given below for these sum rules are valid for 
a general Skyrme EDF given in Eg.® in which all the 
coupling constants could be considered as independent of 
each other. 



After some manipulation we can express them as an inte- 
gral of the dynamical structure function S*'"-* (q, uj) shown 
in Eq. (O as 



(8) 



1. Energy weighted sum rule 

Making the appropriate asymptotic expansion of the 
response functions written in Appendix [BJ we obtain for 
each channel 



Moreover, because of its intrinsic analytic properties, the 
response function satisfies a dispersion relation. As a con- 
sequence, moments can be obtained analytically through 
appropriate expansions in power scries of uj [22| 



for 1 



-oo, the positive odd order moments read 
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x("Hw,q)« 2p^H-(2f+2)A/W^(q), (9) 

p=0 



M. 



(i,±i,i) 



2to* 



1 - 



(1,1) 



4Bt 



(13) 



and can be used for the calculation of the Mi En- 
ergy Weighted Sum Rule (EWSR) and the A/g Cu- 
bic Energy Weighted Sum Rule (CEWSR). 

• for w ^ 0, the negative odd order moments can be 
extracted as 
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X(")(a.,q) « -2pE iuf^M%^,^{q), (10) 

which will be used for the M_i Inverse Energy 
Weighted Sum Rule (lEWSR). 

In the above formula p represents the density of the sys- 
tem. 

The situation is the following : we have at our disposal 
two expressions for M^"' : the first one (Eq. ([8])) is purely 
numerical and implies the whole response function; the 
second one Eq. ^ (or Eq. (jTU]) '). which originates from 
the dispersion relation satisfied by the response function 



If one now takes into account the expression of the 
isoscalar effective mass, i.e. m/m* = 1 -|- 2mpCg, one 
can rewrite the expressions of the EWSR in terms of the 
coupling constants of the Skyrme EDF as 



(S,M,I) 



2m 



(14) 



where a'-^''^'^-* is a sum of contributions corresponding to 
each part of the Skyrme EDF which is considered, cen- 
tral and tensor parts of the Table I for instance. The free 

2 

part ^ corresponds to the kinetic part of the hamil- 
tonian since only the gradient terms of the interaction 
contribute to the corresponding a^^'^'^-' coefficient. Note 
that the spin-orbit part of the Skyrme interaction does 
not contribute to the EWSR. 

In self-consistent RPA calculations, i.e. when the same 
effective interaction generates the HF mean field and 
also produces the residual interaction, positive odd-order 
RPA sum rules can be calculated through the Thoulcss 
theorem, by taking the expectation values of appropriate 



operators on the HF ground state (see for example Bo- 2. Cubic energy weighted sum rule 

higas et al. (30l | and Lipparini et at. [20l | for the details of 
this technique). For the Mf'^'^^ EWSR one can write 



^(S.MJ) ^1^01 
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Q 



(S,MJ) 



|0), (15) 



calculated for each (a) channel with the operator 

Q(S,MJ) ^ giq.r, 0(S,M,I) .^^ ^^^jg J ^.^j^ 

the hamiltonian H built up with the zero range Skyrme 
effective interaction. We have checked, after some te- 
dious calculations, that this result coincides exactly, as 
it should be, with that obtained with the asymptotic ex- 
pansion of the response. Note that the double commu- 
tator technique (Eq. [15]) uses the full hamiltonian H of 
the system with a Skyrme interaction and it can not be 
used with a generalized EDF which does not derive from 
a Skyrme interaction. 



TABLE I: Operators used in each (S, M, I) channel. Columns 
3 and 4 give the central and tensor contributions to the 
EWSR respectively. a'^,af^ {ji respectively) are the stan- 
dard components of the vector CTi defined as erf = af and 
a±i = T('xf ±^af)/V2. 
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Making the expansions of the responses (see 
pcndix|B]), we obtain successively for each channel 
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with the VFj- ' coefficients given in Appendix [Cl and the 
usual relations : p = 2kp/{37T^),k = 

In principle the M^^'^'^^ CEWSR can be also obtained 
analytically from the commutator machinery briefly de- 
scribed in the previous paragraph but the operator to be 



considered involves now a triple commutator, i.e. three 
times the interaction. The calculation becomes then very 
tedious and has been checked only for the central part of 
the Skyrme interaction [3l| which gives a CEWSR which 
does not depend on the value of spin projection M 
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3. Inverse energy weighted sum rule 

This moment cannot be obtained through the commu- 
tator machinery but only using appropriate constrained 



Hartree-Fock calculations through the well-known dielec- 
tric theorem [30j . 

Thus, we use again the appropriate expansion of the 
expressions given in Appendix [B] to obtain finally 
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with /(fc) 
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Since the in- 



2k(i-^^)iog(s); 

stabilities we arc looking for, arc related with poles of the 
response functions at zero energies, this sum rule will be 
shown in the next part to be of fundamental importance 
for the detection and therefore prediction of instabilities. 
Since in the definition of the contribution of the 

low-energy part is more important because of the factor 
l/w in the integrand, this sum rule is more sensitive to 
the poles than the others. 



center of a nucleus effectively explores, because of fluc- 
tuations, not only the saturation density but also some 
larger values for which one may observe a divergence of 
the response functions. In the following, we will concen- 
trate ourselves on the detection of such poles. 



B. Sum rules 

1. EWSR 



III. RESULTS 
A. Response functions 

We have already discussed in article I the fact that 
the tensor may contribute significantly to the response 
functions. Here, we precise several related aspects in view 
of the forthcoming discussion about instabilities. 

Quite generally, in S" = channel, the tensor terms do 
not affect qualitatively the response; all tests performed 
using TIJ tensor interactions discussed in ref. |32| exhibit 
the same qualitative behavior. The situation is quite dif- 
ferent in S" = 1 channels; the effect from the tensor terms 
is large whatever the value of the spin projection M is. 
Actually, depending on the values of the transferred mo- 
mentum q and the density p, the response functions can 
even increase significantly and diverge for finite q for a 
certain critical density pc- As illustrated in Fig. 1, one 
can typically observe two types of extremes phenomena : 
the first one (left panel) corresponds to an accumulation 
of strength at finite energy (and low transfer momen- 
tum) and is related to the zero sound mode whereas the 
second one (right panel) is associated to a pole at zero 
energy (and finite momentum). Although a one-to-one 
correspondence between infinite matter and nuclei is ob- 
viously not trivial, preliminary tests seem to show that 
the latter divergence actually reveals the presence of in- 
stabilities observed in nuclei [iTj , with the appearance of 
domains with typical size of the order of 21: /q (ssj . The 



As an example. Figure [2] shows the EWSR calculated 
for the equilibrium density, for the T44 tensor parameter- 
ization and for all the six (S, M, I) channels. In each case 
the result obtained with the integral (Eq. ^) is com- 
pared to the exact calculation (Eqs. (fTT|) . (fT2|) and p^ ). 
As expected, both results coincide, satisfying then the 
sum rule. It remains an exception for the two (1,0,1) 
and (1, 1,0) channels where the integral calculation vio- 
lates the sum rule. This is actually due to the presence of 
a pole (indicated by an arrow on Fig. 2) in the strength 
function at (/ ~ 1.5 fm~^ for the (1, 1, 0) channel (see cap- 
tion of Fig. 1 right panel) and at g ~ 2.2 fm^^ for the 
(1,0,1) channel. These poles which are clearly exhibited 
in the lEWSR (see below) make the sum rules unphysical 
from and above the q value of the pole. 



2. CEWSR 

For the same example and for the same conditions. Fig- 
ure [3] shows the CEWSR. As for the EWSR we observe 
a perfect correspondence between the two calculations 
of the sum rule, integral or analytical expression (see 
Eqs. (|17p . (fT8|) and p9|) ) except in the channels which 
exhibit a pole in the strength function. The same re- 
marks as for the EWSR can be done. Due to the cubic 
energy weighted in this sum rule, the violation of the 
concerned sum rules does not appear very clearly on the 
figures. 
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(0 [MeV] CO [MeV] 



FIG. 1: Nuclear response function for two given values of the transferred momentum q (in fm"'^) for the T44 tensor parame- 
terization. On the left we show the response function in the channel (1, 1, 1) with the discrete ph transition corresponding to 
the zero sound. On the right we show the response function in the channel (1,1,0) in proximity of a pole.psat is the saturation 
density of the system. 




FIG. 2: (Color on line) EWSR (in MeV) as a function of 
the transferred momentum q (in fm~^) for the T44 tensor 
parameterization. Red dashed (black solid) line correspond to 
the integral (Eq. ((S|) (analytical expressions, Eqs. (|ll|l . p2p 
and (|13|) respectively). Results are shown for the saturation 
density and for each (S, M, I) channel. 



3. lEWSR 

Always for the same example and for the same condi- 
tions, Figured shows the lEWSR. As for the EWSR we 
observe a good correspondence between the integral and 
the analytical expression (see Eqs. (HH), (HH) and 1^^) 
except when there is a pole in the strength function. 
Moreover the discrepancy in that case is very sharp so 



FIG. 3: (Color on line) CEWSR (in MeV^) as function of 
the transferred momentum q (in fm~^) for the T44 tensor 
parameterization. See Figure [2] for other details. 



that the signature of the pole is very clear. This is due 
to the fact that for lEWSR this pole is also present on the 
denominator of the analytical expressions (see Eqs. (|2ip . 
([22]) and (f23|) ). Thus, we show here that there is unique 
correspondence with the pole observed in the strength 
function and the pole of the lEWSR. An immediate con- 
sequence is that Eqs. ((2T|) . (|22p and ([23]) can be used in 
a fit protocol in order to test directly the occurrence of 
instabilities. 

It should also be noticed (see the inset of FigH]) that a 
small amount of strength is missing at low q in channels 
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(1,0,0) and (1, 1, 1). This again corresponds to the zero 
sound mode aheady shown in left panel of Fig. [1] and 
will not be discussed here. 



meet at the critical point. This fact is due to the C^'' 
coupling constant and it can be viewed as a surface effect. 
Without this term in the functional the two branches of 
the spinodal curve would turn into two parallel lines [si 
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FIG. 4: (Color on line) lEWSR (in MeV"^) as function of 
the transferred momentum q (in fm~^) for the T44 tensor 
parameterization. See Figure [2] for other details. 



C. Instabilities 

When the response function exhibits a pole at zero 
energy in a given channel one can suspect that an un- 
physical instability will occur in finite nuclei if the corre- 
sponding critical density is closed to the saturation den- 
sity. The goal of this part is thus to show for typical 
Skyrme parametrisations whether they lead to such prob- 
lem, that is the appearance of an unphysical instability 
when the critical density calculated using the lEWSR 
is close to the saturation density . Since we have shown 
that there is a direct connection between the pole (when 
it docs exist) observed in the response function and the 
pole observed in the M_i sum rule. It is easy to plot the 
critical densities pc as a function of g by simply solving 
l/M_i{pc,q) = in each channel. For example. Fig- 
ure [S] shows the behavior of the critical density for each 
(S, M, I) channel for the interaction T44. As a guide-eye, 
the saturation density Psat is also plotted. As claimed, 
one can clearly see that one exactly obtains the same re- 
sults if one considers the pole of the M^i sum rule (open 
circles) or the pole of the corresponding RPA responses 
(dashed lines) . For this particular paramctrisation, insta- 
bilities appear both in the (1,0,0) and (1,1,0) channels 
at Pc = Psat- For the (0,0) channel, we can also see on 
Figure [5] the well-known spinodal instability. This spin- 
odal instability is viewed here as a two branches curve 
corresponding to the two critical densities observed in a 
standard plot of this spinodal curve. These two branches 



FIG. 5: (Color on line) Critical densities (in fm~"^) as func- 
tions of the transferred momentum q (in fm~^) for the T44 
tensor parameterization. Open circles show the critical densi- 
ties extracted from the pole of the AI-i moment while the 
dashed lines correspond to the pole of the corresponding 
strength function. 

Figure [6] displays the critical densities for some usual 
Skyrme EDF. All the Skyrme EDF exhibit the same 
physical spinodal instability but the behaviors of the crit- 
ical densities in the other channels are very different and 
depend strongly of the parameterization under consider- 
ation. 

Similarly Figure [7] shows the evolution of critical densi- 
ties following two series of parameterizations with tensor 
couplings ni — T16 and Til — T61 when one considers 
the Cq — Cl plane of tensor coupling constants studied 
by Lcsinski et al. [s^ . In this case we only show the 
S = 1 channel for the different TIJ forces, showing that 
they all are unstable against spontaneous polarization of 
finite-size domains as already observed by Hellelmans et 
al. [11. 



D. Landau parameters 

Another important constraint concerning stability of 
a paramctrisation is given by the Landau parameters : 
since they represent the short-range of the interaction, 
they have to be positive. The Fi , spin- independent 
Landau parameters must obey to the stability condition 



Fi 



21 + 1 



> 0. 



(23) 



10 




0.5 
0.4 
0.3 



1 — 

N I M=0 1=0 
S=l M=OI=0 
S=l M=0I=1 
S=l M=l 1=0 
S=l M=l 1=1 




q [fm 



0.5 
0.4 
0.3 
clO.2 



SLy4 



S=0 M=0 1=0 
S=0 M=0I=1 
S=l M=0I=0 
S=l M=l 1=0 





q [fm 



0.5 
0.4 
0.3 



clO.2 



S=0 M=0 1=0 
S=l M=0 1=0 
S=l M=0 1=1 
S=l M=l 1=0 
S=l M=l 1=1 




q [fm 




q [fm" 



FIG. 6; (Color on line) Critical densities, pc(in fm as functions of the transferred momentum q (in fm ^) for some usual 
Skyrme EDF: SkP [sll, SkM* [sl, SGII [l5[, SLy4 HqKH, BSk8 [H, and SkO [IJ. The horizontal dashed-dotted line 
represents the saturation density of the system. 
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FIG. 7: (Color on line) Critical densities, pc(in fm~^), in the S — \ channels as functions of the transferred momentum q (in 
fm~^) for the Til to T16 (top four panels) and for the Til to T61 (bottom four panel) tensor parameterizations.The horizontal 
dashed-dotted line represents the saturation density of the system. 
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Remember that the / = Landau paramctrization can 
be related to the second derivative of the EDF with re- 
spect to 1, 2r , ^cr , 2^(TT ( see [3l for details), the pertinent 
variables of each {S, I) channel. Eq. represents thus 
the condition that the concavity of the equation of state 
(EoS) at the equilibrium must be positive in each {S, I) 
channel. The result is represented on Figure [5] for some 
Skyrme interactions. Similarly the / = 1 Landau param- 
eters can be related to the effective mass in each {S, I) 
channel and Eq. ([25)) requires that each effective mass 
is positive without any pole. Similar conditions exist 
for Gi,G[ spin-dependent Landau parameters. They are 
shown on Figure [9l In the presence of a tensor interaction 
a new additional condition that prevent the deformation 
of the Fermi sphere has to be satisfied. Following the 
derivation of Brown et al. |29| we have 



1 + 3G1 - —Ho > (24) 
1 + + ^Ho > (25) 



1 -f -Gi - -Ho > (26) 
(l + ^)+^V'G§ + 8i/2>0 (27) 

{l + ^)-yGl + 8H^o>0 (28) 

and similarly for Hq. On Figurc[TO]we show the left hand 
side of Eqs. P4II28I1 . This result is consistent with the 
results presented by Cao et al. [TB] , but generalized here 
for the case of a Skyrme functional. 



IV. SUMMARY AND CONCLUSIONS 

In this article, we have presented the analytic con- 
tribution arising from the tensor terms to the RPA re- 
sponse functions with a general EDF as a starting point. 
From these response functions, we derived the Landau 
parameters and we focused on instabilities at zero en- 
ergy and finite transfer momentum. In particular, we 
have shown that a divergence of the response structure 
functions x'"^ (Qi indicates a finite size instability in 
infinite matter. Moreover this instability can be detected 
by simply using the analytical lEWSR, which is a great 
advantage for future applications. At this point, one also 
has to note that a systematic study of the critical densi- 
ties is in progress in order to determine whether the link 



between the divergences of x*^"^ (q, w) and the instabilities 
encountered in nuclei at the Hartree-Fock approximation 
is robust. 

Another important point under study is the identifica- 
tion, directly from the Skyrme energy functional, of the 
contribution of each term of the EDF in the response 
functions. Such a study would enable us to put some 
constraints on the different constants in order to avoid 
instabilities. In the same spirit, a detailed study of sum 
rules can enlighten the contribution of the tensor for var- 
ious physical situations (see for instance [l^]). Finally, 
applications to pure neutron matter is of great impor- 
tance (see for instance [isl - fSTj ) and will be the subject 
of a forthcoming article in preparation. In that case, the 
above formulae are no longer directly usable and have 
been adapted to that specific case. 
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FIG. 8: ^ = Landau parameters as a function of tlie ratio p/ psat, where p is the density of the system and paat is the saturation 
density, for some usual Skyrme parameterizations. 
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FIG. 10: On the left panel we show the left hand side of Eqs. (|24l28p . while on the right panel, adopting the same color code, 
the same but replacing in Eqs. (|24I28I) with the terms Go, G'l, Hq. See text for details. 
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Appendix A: Particle-hole matrix elements of the 
zero range tensor part of the interaction. 



tion for the tensor part of the functionaL 



Following the notation adopted in article I, we give in 
Table the values of the particle-hole residual interac- 



TABLE II: Contribution of the EDF tensor part to the residual interaction in terms of the Bj coupling constants. For the sake 
of simplicity we have introduced the notation: Kij = [(fci2)i(fci2)j], where (fci2)^^' is defined in Eq.(9) of article I. The term 
Sss'^siSii'Sqqi is implicit everywhere. 





M' = 1 


M' = 




M' = -1 


M = 1 


-2g^ {Bf + 4Bf 
+4Bf Ko.o 
-4 {2Bj + Bf ) Ki,_i 


-4Bf K-1,0 




-4B[ K-1,-1 


M = 


4Bf Ko,i 


-2 [Bf - ABf" + 4Bf ' 
+4 (Bf + Bf) E 
-8Bj Ki,_i 


+ B[)q' 

'■0,0 


4 Bf 


M = -1 


-4B[ Ki,i 


-4Bf Ki,o 




~2q^ {Bj + 4B^') 
+4 Bj Ko.o 
-4 {2Bj + Bf) ]Ki,_i 



r 



Appendix B: RPA nuclear responses 

We recall here the nuclear responses already given in 
article I but rewritten here in terms of the coupling con- 



stants of the Skyrme EDF. We keep in mind from article 
I the definitions of these coupling constants in terms of 
the parameters of the Skyrme interaction. 



J 



• For the S = Q channel 



(0,1) ^ 1 

Xrpa 



Xo 



(0,1) 



r;42 74 , *2 22 12* 1 

[kpX2 - kpXoXi + m u} Xo - jl PXo\ 



(0,1) 



1 - \m*pW!^ 



(0,1) 



Xo 



(Bl) 



with 



~(o,i) 



PFf + 16 [Cj 



/S2-/33 



1 + g2 _ /33) + 4Bf - 2Bf 



(B2) 
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• and for the S* = 1 channels 



XHF 
(1,0.1) 

Xbpa 



[l + m*pBf] + 



;(i,o,i) 



Xo 



(1,1) 



2X0 



(1,1) 



4B 
4B 



[1 + 2m*pB[] xo - 24X2 + 2m* pklBf (xo - X2)} 

2 

I 1 i 2 ;4 I *2 22 12* 1 



W^^''^ + ABj + 4B[ 


[1 + imVX(i'O'i)] 


1 — ^m* p 


W!^^'^^ + ABj + ABf - X(i.o.i) 



(B3) 



XHF 

(1,±1,I) 
Xrpa 



= [l-\m*pBf]-W, 



;(1,±1,I) 



Xo 



(1,1) 



w. 



2X0 



(1,1) 



ABJ + 2Bf 



4Bf + 2Bf 



{\q^ [1 - m*pBf] xo - 24x2 - m*pklBf (xo - X2)} 
2 

r;42 74 I *2 22 1*2 1 

{kpX2 - kpXoXA + Xo - P9 Xoj 



(1,1) 



ABj 



[1 



±1,1)1 



1 - \m*p W!^^'"^^ + ABj - iX(i^±i.i) 



(B4) 



were we have used 



;(1,0,I) 



2 
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2q^Bj 
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m u! 



m 



■p[Bf\ 



(B5) 



;(1,±1,I) 



q m p - 



Sq' [C] 



m uj 



VJl 



(fi2 - /33) 



1 + g2 - /33) w, 



Xo 



_(0,I) 

q 



Bi 



X2 +4fc^X4 



(B6) 



The X^-^'^'^") coefficients occurring in the previous ex- Appendix C: the Wj;°'\ W^"^ and X'"' coefficients, 
pressions of the 5 = 1 response functions are defined 
in the Appendix [Cl while the momenta were already 
defined in Appendix D of article I. 

In order to simplify all the written formula in the pres- 
ence of a tensor part in the Skyrme interaction, the w'f''^^ 
and ^2^'^^ coefficients have been defined as 
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and the same way for the coefficients can be 

written in terms of the EDF couphng constants as 



9^ 



-w. 



(1:0) _ 



(14) 



4^ 
^0 ' 

^1 : 



for the residual interaction of a general Skyrme functional 
in the Landau limit (see Eq.([5|)) we define the W^^j^^ and 



^2 L^'' coefficients as 



(14) 

1, L 

(0,0) 

2, L 



2Co^" + (2 + 7)(l+7)Co^V, 
2Cf° + 2Cf'>T, 
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